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I.  INTRODUCTION 


One  measure  of  the  “goodness"  of  an  antenna  or  sensor  array  is  the  sidelobe  structure 
of  the  beampattern  function  when  the  conventional  linear  beamformer  is  used.  In  this  report, 
we  examine  the  performance  of  the  adaptive  beamformer  in  terms  of  spurious  peaks  in  the 
array  response  function  by  exploiting  the  relation  between  adaptive  array  response  and  the 
conventional  beampattern.  An  expression  that  provides  an  approximate  relationship  between 
the  array  response  function  for  an  adaptive  beamformer  and  the  conventional  beampattern 
function  in  the  sidelobe  region  is  developed.  The  validity  of  the  approximation  is  examined 
using  simulations  for  a  variety  of  arrays  including  deterministic,  randomized,  and  random 
cases.  The  approximate  relationship  is  used  to  obtain  the  probability  that  a  spurious  peak  in 
the  adaptive  array  response  exceeds  a  threshold  for  random  arrays. 

The  organization  of  the  report  is  as  follows:  In  Sections  11  and  III,  the  basic  relation¬ 
ships  for  conventional  and  adaptive  beamforming  are  established.  In  Section  IV,  the  approxi¬ 
mate  relationship  between  adaptive  array  response  and  conventional  beampatterns  is  deve¬ 
loped.  In  Sections  V  and  VI,  the  relationship  is  validated  through  simulation.  In  Section  VII, 
the  spurious  peak  estimator  for  random  arrays  is  developed. 
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II.  CONVENTIONAL  BEAMFORMING  AND  BEAMPATTERNS 

The  relationship  between  the  beampatterns  or  power  patterns  of  an  array  and  the 
observed  array  response  patterns  when  using  conventional  beamforming  is  fairly  straightfor¬ 
ward.  For  example,  when  the  noise  field  consists  of  a  single  signal  from  a  (multidimensional) 
arrival  angle  6,  the  conventional  array-response  pattern  (normalized  to  0  dB  at  0)  would  simply 
be  the  conventional  beampattern  steered  at  the  signal 

I^J^lookl^  • 


Here,  the  vector  S  represents  an  ideal  plane  wave  and  has  the  form 


S 


1 


exp 


OJ 

—  <  Z|,  k:  >) 

c 

OJ 

-  <  Zw,  K  > 

C  ‘ 


(2) 


The  convention  in  this  report  will  be  to  use  unit  norm  signal  vectors.  Frequently,  norm-N 
vectors  are  used.  When  necessary,  a  factor  of  N  will  be  explicitly  included.  The  unit  vector  k 
indicates  the  direction  of  look  or  steering  angle  as  well  as  the  direction  of  arrival  of  the  signal. 
Element  position  vectors  for  the  N  elements  of  the  array  are  Z]  through  Zj^.  The  inner  product 
is  represented  by  <•,•>.  Vectors  and  matrices  are  denoted  by  bold-face  characters  and  • 
denotes  complex-conjugate  transpose. 


More  generally,  we  assume  that  the  time-series  data  at  the  array  can  be  written  as  a 


vector  X(t)  = 


Xi(t) 


.  Let  X(a>)  represent  a  Fourier  transform  vector 


T 

X(w)  =  —  f  X(t)e-‘"*  dt  , 
T 

o 


(3) 


where  the  frequency  a>  is  an  integral  multiple  of  2}r/T.  The  data  at  frequency  w  are  then 
completely  described  by  a  cross-spectral  matrix  R  which  is  formed  by  summing  outer  products 
of  Fourier  transform  vectors 


R((u) 


—  X  Xg(a,)x*g(a»)  . 
^  8=1 


(4) 


If  the  signal  plus  noise  field  consists  of  K  uncorrelated  signals  Sj  of  power 
o?  (i  =  1,2,...,K)  in  noise  of  power  o?  which  is  uncorrelated  from  sensor  to  sensor,  then 

S:  ' 


2 


(5) 


K 

R  =  S  Na2  SjS*  +  oj, 
i=l 

Another  way  to  write  this  is 

R  =  ADA*  +  ojl  ,  (6) 


where  A  is  the  N  x  K  matrix  of  discrete  signals 


and  D  is  the  diagonal  matrix  of  signal  powers 


D  = 


No: 


No: 


n4 


The  conventional  beamformer  calculates  S*(jo|jRS|qqJj  and  peaks  in  the  resulting  array- 
response  pattern  are  identified  as  signal  arrivals.  The  power  output  of  the  conventional  beam- 
former  with  the  signal  plus  noise  field  of  (6)  is 


^CBF  = 


S*ook®S 


look 


=  S;„„kADA*S,^K  »  oj 

i®Iookl^ 


=  1*  Sfljisjs 

L  i=l 


(7) 


with  /3j  =  No^./o^-  resulting  conventional  array-^response  pattern  is,  thus,  heavily  depen¬ 
dent  on  the  conventional  beampatterns  steered  at  the  signals.  Figure  1  is  the  17-clement  array. 
Figure  2  is  the  conventional  beampattem  for  the  I7-elemcnt  array  of  figure  1  where  it  is  steered 
to  9  =  180®  in  uncorrelated  noise.  Figure  5  is  the  conventional  array  response  which  results 
when  a  single  0-dB  signal-to-noise  ration  (SNR)  signal  is  present  in  uncorrelated  noise.  Note 
that  the  array-response  plot  is  a  scaled  version  of  the  beampattem.  The  sidelobes  in  the  beam- 
pattern  become  spurious  peaks  in  the  array-response  plot. 
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III.  ADAPTIVE  BEAMFORMING 

Adaptive  beamformers  rely  on  information  about  the  noise  field  to  imp:  ove  signal 
detection  performance  and  arc  capable  of  extremely  high  resolution,  low  bias,  and  interference 
rejection.  In  some  cases,  adaptive  beamformers  lack  the  robustness  of  conventional  beam- 
formers.  Nevertheless,  the  strengths  of  adaptive  techniques  are  considerable,  and  they  offer  a 
powerful  tool  for  target  detection  and  localization. 

One  outstanding  issue  is  the  impact  of  an  array’s  conventional  beampatterns  on  its 
adaptive  array  response.  One  might  argue,  on  purely  heuristic  grounds,  that  an  array  that 
performs  poorly  using  conventional  techniques  will  also  have  relatively  poor  adaptive  perfor¬ 
mance.  In  fact,  this  is  the  case.  However,  we  wish  to  make  this  relationship  more  quantitative 
and  relate  it  to  signal  detection. 

We  consider  two  adaptive  beamformers.  The  first  is  the  minimum-energy  (Capon’s) 
beamformer.  This  is  closely  allied  with  the  maximum-entropy  beamformer  described  by  Burg. 
The  second  is  the  music  or  eigenvector  beamformer  which  is  currently  a  topic  of  much  interest. 

Given  a  noise  field  described  by  a  cross-spectral  matrix  R,  the  minimum  energy  beam- 
former  selects  a  filter  vcct<^r  which  minimizes  the  output  power  subject 

to  the  constraint  that  W  =  1  where  is  an  ideal  plane  wave  in  the  look  direction. 

The  solution  to  this  constrained  optimization  problem  is  found  by  the  method  of  Lagrange 
multipliers.  We  write 

F  =  W*RW  +  A(W*S  -  1)  . 

and  set 

dF  dF  3F 

-  =  0  , -  =  0  ,  and  —  =  0  . 

aw  ax 

The  familiar  minimum  energy  steering  vector  is  thus 

R~‘s 

W  =  — -  ,  (8) 

S  R~*S 


and  the  output  power  of  the  minimum  energy  beamformer  is 


^ME 


W*RW  = 


"r-‘s^ 

•  R  . 

"r-‘s'" 

1 

[s'r->sJ 

®look*^  ^®look 

(9) 


Return  to  the  noise  model  of  (6),  R  =  ADA*  +  o^I,  and  assume  that  the  number  of 
signals  K,  is  less  than  the  number  of  arr^  elements  N.  We  note  that,  provided  some  basic  array 
design  criteria  is  met,  both  A  and  ADA  are  rank  K  matrices.  Also,  since  D  is  positive-definite, 
ADA  is  positive  semi-definite.  We  can  then  write  ADA  in  terms  of  its  eigenvalues  and 
eigenvectors  as 


4 


(10) 


K 

ADA*  =  X  A.jVjV* 
i=l 

where  the  eigenvalues  A.,  are  all  strictly  positive.  The  space  spanned  by  the  K  eigenvectors  of 
ADA*  associated  with  positive  eigenvalues  is  also  the  space  spanned  by  the  columns  of  ADA* 
—  the  signal  vectors.  This  space  is  called  the  signal  subspace. 

With  the  result  of  (10)  and  the  fact  that 
N 

I  =  5;  ViV’ 

i=l 

for  any  set  of  orthonormal  vectors  {Vj}  of  length  N,  we  can  write 
K  N 

R  =  2  *  "Sivl  *  s  (M) 

i=l  i=K+l 

Thus,  the  K  largest  eigenvectors  (those  associated  with  the  K  largest  eigenvalues)  form  a  basis 
for  the  space  spanned  by  the  K  signal  vectors.  The  (N-K.)  smallest  eigenvectors  are  orthogonal 
to  the  signal  subspace  and  form  a  basis  for  the  noise-only  or  noise  subspace.  Furthermore,  the 
(N-K)  smallest  eigenvalues  are  all  equal  to  the  noise  power  o^. 

The  second  of  the  two  adaptive  beamforms  we  consider,  the  music  beamformer 
(Johnson.  September  1982,  and  Schmidt,  1979)  takes  advantage  of  the  orthogonality  of  the 
noise-only  eigenvectors  with  the  signal  vectors.  The  output  of  the  music  beamformer  has  the 
form 


^MUSIC  = 


'look 


N 


2  -ViV* 


vfK+1 


*look 


(12) 


Peaks  in  the  MUSIC  response  plot  corresponding  to  signal  directions  of  arrival  occur  when 
Siook  Is  in  the  signal  subspace  and  the  denominator  of  (12)  is  zero.  Equation  (12)  is  Johnson’s 
(August,  1982)  implementation  of  this  beamformer. 

Before  investigating  the  array  response  of  these  adaptive  beamformers,  we  look  at  the 
adaptive  beampattems.  An  adaptive  beampattem  results  when  a  filter  vector  Wg^g^r  such  as 
the  vector  of  (8)  is  fixed  and  I  is  plotted  as  we  sweep  an  ideal  plane- wave  S  through 

all  directions  of  arrival.  The  beampattem  is  a  visual  picture  of  the  relative  sensitivity  of  the 
array  to  various  arrival  directions  while  steered  with  Vf steer 

For  the  conventional  beamformer,  the  beampattem  is  independent  of  the 

noise  field.  Since  the  selection  of  the  adaptive  steering  vector  depends  on  the  noise  field, 
however,  this  is  not  the  case  for  the  adaptive  beamformer.  The  beampattem  will  be  heavily 
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dependent  on  the  noise  field.  Figure  4  is  the  beampattern  of  the  minimum-energy  (ME)  beam- 
former  of  the  array  of  figure  1  steered  to  =  180°  when  a  0-dB  SNR  signal  is  present  at  =  30° 
in  uncorrelated  noise. 

At  least  one  property  of  adaptive  beampatterns  can  be  explored,  however.  Return  to 
expression  (8y  for  the  minimum-energy  steering  vector.  If  we  write  the  spectral  decomposition 
of  R  as 


N 


R  =  X  ^iViV  •  , 
i=l 

with  X[  >  X  >  .  .  .  >  and  recall  that 


N 


"■'  =  S  -  ViV'  , 


i=l  ^i 


then  the  beampattern  becomes; 


l<eerS|^  =  with  m  =  - 


®  steer®  ^^steer 


= 


r 


r 


V. 


U 

i=l  '^1 


•N 


VV* 

1 


* 

G 

*^steer 

J 

S 

N 


^  “i  • 

2  —  VjS 


i=l  ^ 


with  Oj  - 


Note  that  for  the  case  of  K  signals,  Xj  »  i=l,2....,K.  Thus,  the  contribution  to  the  beam- 

pattern  due  to  the  signal  subspace  eigenvectors  is  much  less  than  that  due  to  the  noise 
subspace.  We  can  then  write  that 


|W*  rS|2  - 


N  • 

^  «i  • 

i=K+l  ^ 


If  we  look  in  the  direction  of  one  of  the  signals  (i.e.,  let  S  =  S^jg^gj),  then  the  adaptive  beam- 
pattern  will  have  a  value 


2  ^  „2 


l^steer^signall  ~  ^ 


N 


i=K+l 


signal 


=  0 


due  to  the  orthogonality  of  the  signals  with  the  noise  subspace. 


(13) 
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This  is  the  important  principle  of  null  placement.  The  adaptive  beamformer  tries  to 
place  nulls  in  the  beampattern  in  the  direction  of  the  interfering  signals.  This,  of  course,  is  how 
the  minimum-energy  beamformer  minimizes  output  power.  Note  in  figure  4  that  the  ME 
beampattern  has  a  value  of  — 66  dB  at  d  =  30°.  In  fact,  if  the  signals  are  of  infinitely  high 
SNR,  the  nulls  will  be  perfect.  The  MUSIC  beamformer,  by  excluding  the  signal  space  eigen¬ 
vectors  from  its  projection  matrix 

^  I 

2  -  Viv;  • 

i=K+l  ^ 


effectively  increases  the  SNR  of  the  signals  and  the  nulls  are  perfect  regardless. 


IV.  ADAPTIVE  ARRAY  RESPONSE  AND  APPROXIMATION 

Next  we  will  investigate  the  array  response  of  the  adaptive  beamformers.  Once  again 
return  to  expression  (8)  for  the  minimum  energy  steering  vector.  With  n  and  the  aj’s  as  before, 
we  have 


Steer 


N 


S  -ViV* 


i=l  ^i 


’steer 


N 


N 


2  — 


i=l 


a\- 
“r  1 


N 


'N 


N 


S  —  «iV,  -  X  «iVi  ^ 

i=l 


i=i  H 


steer 


’steer 


-  2 


“iVi 


i=l 


(14) 


If  the  signal  plus  noise  field  is  that  described  by  (6),  then  we  can  take  advantage  of  the 
multiplicity  of  the  smallest  eigenvalue  (A.j^+|  =  Xj^+2  ”  •  •  •  “  “  ®n) write 


W 


steer 


H-oi 

^  I  il 


’steer 


i=I 


«iVi 


(15) 


Also,  the  K  largest  eigenvectors  span  the  signal  space.  With  Sj  representing  the  i*^  unit  norm 
plane-wave  signal  (i  =  1,2,  ...,  K),  we  can  perform  a  change  of  basis  and  write 


r  ^ 

^steer  ^steer  “  2 

L  i=l  , 


(16) 


where  the  Wj’s  are  (as  yet  undetermined)  complex  weights. 

We  proceed  to  derive  an  approximation  for  the  adaptive  array  response  in  a  field  of  K 
discrete  signals  and  uncorrelated  noise  of  power  a^.  Since  W^^g^j.  =  ^R'^Sj^ggp  we  see  that 

®  *^steer  ^ 

4 


’steer 


K 

2 

i=l 


WjSj 
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The  output  power  of  the  minim unvenergy  (ME)  beamformer  is  then  ^ven  by 


ME 


<1* 

^  Steer  steer 


1 


’steer 


K  ^  ^ 

steer  '  S  ^^i^i 
V-  i=l 


1=1 


The  output  power  of  the  MUSIC  beamformer  is  given  by 

I 


MUSIC  - 


r  N 


’steer 


X  -  VjV' 
vi=K+l  ^  ^ 


’steer  • 


In  a  similar  fashion,  we  derive: 


Pmusic 


’steer 


v. 


N  7 
°n 


1 

i=K+l 


®steer 


i=l 


J 


withoj  =  VjSstee, 


I 


r 

1 

K 

^steer 

^steer  “  X 
i=l 

J 

4 


'  -  X  Wis;,„,Si 
i=l 

Which  is  the  same  as  equation  (17). 
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z  W 


One  way  to  evaluate  the  weights  W  would  be  to  calculate  the  eigenvertors  and  eigen¬ 
values  of  R  and,  given  the  signal  vector  Sj,  change  bases  from  (15)  to  (16).  This  would  be  cum¬ 
bersome,  however,  and  little  insight  would  be  gained.  Therefore,  we  take  another  tack.  For  any 
given  steering  direction  we  require  that  perfect  nulls  be  placed  in  the  adaptive  beam- 

pattern  in  the  direction  of  each  of  the  K  discrete  signals.  That  is,  we  require  that 


2 

''  K  ^ 

''^steer^j 

= 

oi 

s’sleerSj  ^  1  w’s'Sj 
L=1 

J 

for  j  =  1,2, ...,  K.  This  gives  rise  to  a  system  of  equations; 


I  s;s2  . 

s;sk' 

W, 

®l®steer 

S2S2  1 

W2 

= 

^2^steer 

_  sjs,  . 

.  .  1 

1 

1 _ 

_  ^K^steer  _ 

(18) 


In  another  notation: 


A*A  * 


w. 


=  A*s 


steer 


(19) 


[wkJ 

where  A  is  the  N  x  K  matrix  whose  columns  are  the  signal  vectors,  as  before.  The  weights  are, 
thus,  a  function  of  steering  angle  and  are  related  to  the  conventional  beampatterns. 


With  the  weights  from  (19),  the  minimum  energy  and  MUSIC  response  plots  are 
approximated  by 


^ABF 


o 


2 

n 


1  ®steer'^(^  ^  ^steer 


(21) 


Looking  at  A*A  in  expression  (18),  we  see  that  if  the  signal  vectors  are  mutually  orthogonal, 
the  cross  terms  are  zero  and  we  have  the  result  that  Wj  =  In  general,  they  are  not.  The 

S  jSj  represent  sidelobe  levels  in  the  conventional  beampatterns.  Provided  the  source  arrival 
directions  are  greater  than  a  beamwidth  apart,  however,  the  magnitude  of  the  cross  terms  is 
small.  The  terms  are  small  in  the  same  sense  that  interelement  correlations  in  isotropic  noise 
are  small  if  elements  are  a  correlation  length  apart.  Thus,  we  assume  that  A  A  is  a  diagonally 
dominant  matrix  and  (A*A)~^  =*  I. 
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Finally,  we  have  an  approximation  for  the  adaptive  array  response  (valid  in  the  side- 
lobe  region): 

Pabf  “ - - -  •  (22) 

1  -  X  l^i^steerl^ 
i=l 

This  relates  the  adaptive  array  response  to  the  conventional  beampatterns  steered  at  each  of  the 
discrete  signals. 


V.  VALIDATION  OF  APPROXIMATION 


In  a  series  of  simulations,  expression  (22)  was  validated.  Figures  6  through  19  show 
plots  of  the  actual  ME  array  response  and  the  adaptive  beamforming  (ABF)  approximation  (22) 
for  the  20-element  half-wavelength  line  array  of  figure  5.  Two  signal  plus  noise  fields  were  used. 
In  the  first,  a  single  20-dB  SNR  signal  was  present  at  S  =  30°  in  uncorrelated  noise.  For  the 
second,  three  20-dB  SNR  signals  ^ere  present  at  0  =  30°,  90°,  and  135°.  Included  in  the  plots 
are  some  conventional  beampatterns  steered  at  discrete  signals  normalized  to  the  output  power 
in  the  direction  of  the  signals.  Also  included  are  details  of  the  array-response  pattern  in  the 
sidelobe  region.  Notice  in  figure  7  the  high  level  of  agreement  between  the  actual  ME  array 
response  and  our  ABF  approximation.  Also  note  how  closely  the  ABF  response  resembles  a 
scaled  version  of  the  conventional  beampattern  of  figure  6.  MUSIC  beamformer  array- 
response  plots  are  not  shown  in  this  series,  but  it  was  found  that  MUSIC  array  responses  were 
nearly  identical  to  ME  array  responses  in  the  sidelobe  region. 

Figures  10  and  1 1  are  a  set  of  plots  for  a  20-element  quarter-wavelength  line  array  in 
the  second  noise  field.  Here  we  also  see  good  agreement  between  the  actual  ME  array  response 
and  the  approximation.  Figures  12  through  16  comprise  another  set  with  the  same  noise  fields 
and  the  17-element  Mill's  cross  array  (figure  1), 

Of  particular  interest  is  figure  15.  A  high  spurious  peak  in  the  array-response  plot 
occurs  at  =0  =  270°.  The  approximation  appears  higher  than  the  actual  peak  here.  We  also  see 
in  this  figure  that  the  conventional  beampattern  steered  at  the  second  signal  (0  =  90°)  has  a 
high  sidelobe  at  0  =  270°.  This  sidelobe  manifests  itself  as  a  high  peak  in  the  array-response 
plot  at  that  location.  Thus,  this  peak  is  a  characteristic  of  the  array  and  not  of  the  noise  field. 

The  final  array  type  used  in  this  series  was  a  random  planar  array  in  which  the 
elements  were  distributed  uniformly  on  a  circle.  Figure  17  shows  a  random  circle  of  diameter 
6  wavelengths.  The  conventional  beampatterns  in  the  plots  which  follow  show  a  much  different 
structure  than  those  of  the  line  arrays  or  the  deterministic  Mill’s  cross  array.  The  patterns  are 
much  less  regular  and  the  sidelobes  less  pronounced.  This  is  a  characteristic  of  random  arrays. 
Figures  18  through  21  are  in  a  series  identical  to  those  previous  for  the  random  circle. 

It  is  instructive  to  review  the  assumptions  which  led  to  expression  (22),  an  approxima¬ 
tion  for  the  ABF  array  response  in  the  region  away  from  the  signals.  The  first  was  that  the 
signal  plus  noise  field  consists  of  K  signals  in  uncorrelated  noise  of  power  The  second  was 
that  the  nulls  in  the  beampatterns  are  perfect.  The  third  was  that  (A*A)~*  =  I.  This  last  is  most 
open  to  question.  For  one  thing,  expression  (22)  now  has  a  singularity  wherever 

K 

S  >  1  (generally  near  the  signals). 

i=l 

As  a  bound  on  ||(A*A)~*  -  I||,  we  have 


III  -  (A*A)->||  < 


111  -  A*A|| 

1  -  III  -  A*A|| 


(23) 
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(Appendix  I).  As  a  typical  value  for  |  |I  -  A  A|  1  consider  that  for  a  random  array,  |S  are 
random  variables  with  mean  1/N  and  variance  1/N^  provided  the  signals  are  greater  than  a 
beamwidth  apart.  Then, 


III 


A*  A I 


,  K  K 

-IX 

K  i=l  j=l 


|SiSj|2 


and 


III  -  (A*Ar'ii  < 


For  K  =  3  and  N  =  20,  this  value  is  —  .46. 
We  can  rewrite  expression  (22) 

Pabf  “ - - - 

i=l 

by  writing  it  as  a  geometric  series; 


K 

^  K  ^ 

2 

^ABF  ^ 

1  *  X  l^i^steerl^ 

X  l^i^steerl^ 

+  . . . 

i=I 

i=l 

^  J 

This  could  be  truncated  for  small  to 


PaBF 


4 


1  + 


K 

X  |S>, 

i=l 


steer' 


(24) 


This  is  very  similar  to  expression  (7)  for  the  conventional  beamformer  output  power.  Regard¬ 
less  of  how  we  express  it,  however,  it  is  clear  that  the  adaptive  beamformer,  as  well  as  the 

K 

conventional,  has  a  spatial  spectrum  which  is  dependent  on  the  term  X  l^^^steerl^' 

i=I 
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VI.  AN  EXAMPLE  USING  A  RANDOMIZED  ARRAY 


A  particularly  dramatic  example  of  the  impact  of  ISjSjjgg^r  on  the  adaptive  beam- 
former  is  shown  in  the  following  series  of  plots.  Figure  22  is  of  a  20-element  rectangular  array. 
Figure  23  shows  the  conventional  beampattern  steered  at  6  -  15°.  There  is  an  exceptionally 
high  sidelobe  at  ~d  -  167°.  Figure  24  shows  the  ME  array  response  of  this  array  to  a  signal 
plus  noise  field  consisting  of  a  single  0-dB  SNR  signal  aid  -  15°  in  uncorrelated  noise.  The 
tremendous  (about  7.5  dB)  spurious  peak  at  —6  -  167°  arises  from  the  array.  It  would  certainly 
be  counted  as  a  false  detection  if  not  recognized  as  such. 

As  previously  suggested  in  comments  on  the  random  circular  array,  beampattems  of 
random  arrays  fail  to  display  to  extremely  high  sidelobes  which  arise  due  to  the  symmetries  in 
most  deterministic  arrays.  By  randomizing  an  array  (i.e.,  adding  a  random  component  to  the 
element  position  vectors)  peak  sidelobes  can  be  reduced.  The  effect  on  signal  detection  is  that 
spurious  peak  levels  are  decreased. 

With  this  in  mind,  the  deterministic  positions  Z  = 
figure  22  were  perturbed  by  i  1.  The  x  and  y  random  variables  (measured  in  wavelengths) 

Ly  J 

were  selected  uniformly  on  the  interval  [-t/2,  +t/2].  Selection  of  the  parameter  e  allows  one  to 
move  between  the  deterministic  array  (t  =  0)  and  the  totally  random  array  (t  —  .5  in  this  case). 
Figure  25  through  27  show  the  conventional  beampattems  steered  to  0  =  15°  as  the  randomiza¬ 
tion  parameter  t  is  increased  from  0.1  to  0.5.  For  small  perturbation  (e  =  0.1, 0.3), the  peak  at 
0  -  167°  is  still  present  although  with  somewhat  less  magnitude.  At  higher  values  (c  =  0.4,  . . .), 
the  array  is,  effectively,  totally  random  and  the  deterministic  peak  is  submerged  in  the  random 
pattern. 

Figures  28  through  30  show  the  corresponding  ME  array-response  plots.  For  the  deter¬ 
ministic  and  slightly  randomized  arrays,  the  spurious  peak  is  pronounced  but  decreases  rapidly 
as  the  array  is  more  random.  Figure  31  displays  graphically  the  correspondence  between  peak 
sidelobe  level  (in  the  conventional  beampattern)  and  peak  spurious  peak  level  (in  the  array 
response).  It  also  displays  the  effectiveness  of  using  randomization  to  decrease  both. 


^  yj  of  the  rectangular  array  of 
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VII.  SPURIOUS  PEAK  ESTIMATOR  FOR  RANDOM  ARRAYS 


For  an  array  in  which  the  element  position  vectors  are  random,  the  quantity 


N  N 


2  2 


i  —  (»fj  'fsteer^n  ^m^ 
e 


(25) 


is  a  random  variable  provided  S|  and  are  greater  than  a  beamwidth  apart.  It  has  a 
probability  density  function  given  by 


(26) 


where  =  1/N  (Steinberg,  1976,  eq.  8.18).  This  density  has  mean  1/N  and  variance  l/N^  and 
represents  the  probability  that  the  sidelobe  level  will  attain  a  certain  value. 

If  we  assume  that  the  beampatterns  are  independent  of  each  other,  then  the  sum 


K 

y  =  2  l^i^steerl 

i=l 


2 


has  a  probability  density  function  g(y)  given  by 

g(y)  =  -  yK-1  .  y  >  0  (Appendix  2) ,  (27) 

(K  -  l)!(o2)K 

The  output  power  of  the  adaptive  beamformer  for  a  random  array  is,  thus,  a  random 
variable  in  the  region  away  from  the  signals.  The  transformation  is 


1  2)  ISjSjtggfP 


i=l 

The  probability  that  y  is  greater  than  or  equal  to  )  and  tne  output  power  is  infinite  away  from 
any  of  the  signals  is 
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which  is  negligible  for  most  reasonable  choices  of  K  and  N  >  K.  The  output  power  P in 
the  region  away  from  the  signals,  has  a  probability  density  function  h(p)  given  by 


h(p)  = 


4 


p2(K  -  l)!(o^)‘‘ 
and  a  distribution  function 
P 

h(p)  =  J  h(t)dt 
^2 


— 1 

_ 1 

K-I 

1 

L  "  - 

exp 

P 

L  J 

:p^<^  (28) 


1 

4^ 

1 

4 

=  1  -  exp 

-  — 

1  - 

• 

1  +  — 

1  - 

_  a2 

L  p  J 

L  a2 

L  p  J 

r 


1 


(K  -  I)! 


1 


1  - 


K-l 


_i 


(29) 


We  develop  our  estimate  for  a  peak  spurious  peak  in  much  the  same  way  that  Stein¬ 
berg  (1976)  arrives  at  a  peak  sidelobe  estimator.  Since  the  bcampatterns  of  an  array  are 
completely  described  by  some  number  of  independent  samples  n  (determined  from  the  correla' 
tion  properties  of  the  beampattern,  and  the  array-response  plot  is  a  function  of  those  beam- 
patterns)  it  follows  that  the  array-response  pattern  is  described  by  n  samples  also. 

For  a  random  circular  planar  array  of  diameter  D,  the  correlation  structure  is 
described  by 


where  u  is  related  to  the  angle  between  look  directions  y  by 


u  =  ±2  sin 


and  j|  is  the  first  order  Bessel  function  of  the  first  kind.  The  number  of  independent  samples  or 
array  parameter  n  is  given  by 


t 
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n  =  -  .  (30) 

sin"*  1^.60984 

The  probability  that  a  peak  spurious  peak  will  not  exceed  some  threshold  /3  is  then  given  by 


Pr[peak  spurious  peak  <  )3]  =  [H()3]”  ,  (31) 

where  H(/9)  and  n  are  given  by  (29)  and  (30). 

A  number  of  20-element  random  circular  arrays  of  radius  D  =  lOX  were  generated.  In  a 
noise  field  consisting  of  0-dB  white  noise  and  three  signals  with  various  arrival  and  strength 
structure,  minimum-energy  array-response  plots  were  calculated.  In  50  trials,  a  peak  spurious 
peak  was  found  and  a  cumulative  histogram  of  these  values  appears  on  figure  32.  Also  plotted 
is  [H(^)]"  against  the  threshold  in  dB  for  the  appropriate  parameters. 


rr 

-  =  51 

sin~^[.60984/10] 


1 


1 

20 


K  =  number  of  signals  =  3  . 


(32) 


This  shows  remarkable  agreement! 

Figure  33  displays  a  series  of  [H(/8)]”  against  p  for  a  variety  of  20-clement  random 
circle  arrays.  The  signal  plus  noise  field  is  assumed  to  consist  of  three  signals  in  0-dB  uncorre¬ 
lated  noise.  The  independent  variable  is  the  array  parameter  n  or,  equivalently,  the  array  size. 
For  a  fixed  number  of  elements,  an  increase  in  diameter  while  decreasing  mainlobe  width  (and 
thus  increasing  resolution  capabilities)  also  has  the  effect  of  raising  the  probability  of  spurious 
peak  threshold  exceedance  (thus,  increasing  the  probability  of  false  detections).  Peak-sidelobe 
estimates  for  the  same  arrays  are  displayed  in  figure  34.  The  correspondence  between  peak 
sidelobes  and  peak  spurious  peaks  is  clear. 

Figures  35  through  38  illustrate  the  ME  array-response  plots  for  some  realizations  of 
20-element  circular  random  arrays  of  various  diameters.  The  signal  plus  noise  field  consists  of 
three  signals,  20-dB  SNR  at  «  =  0°,  10-dB  SNR  at  0  =  45°,  and  0-dB  SNR  at  0  =  90°,  in 
uncorrelated  noise.  Figure  35  also  shows  the  MUSIC  array  response  and  the  conventional 
array  response.  Note  the  appearance  of  spurious  peaks  away  from  the  signals  in  the  adaptive 
array  responses. 

For  a  random  line  array  of  length  L,  the  correlation  structure  is  described  by 
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L  ^  j 

w 

where  u  is  related  to  the  steering  angle  d  measured  from  broadside  by 
u  =  sin  0  -  sin 

Here,  0^  is  the  steering  angle.  The  array  parameter  n  is 
L 

n  =  —  [1  +  |sin  0  |] 

A 

The  array  parameter  is,  thus,  dependent  on  steering  direction  (or,  in  the  context  of  signal  detec¬ 
tion,  signal  direction  of  arrival).  Since  the  signals  arrive  from  all  directions,  we  require  a  repre¬ 
sentative  value  for  n.  For  this  we  use 


Figure  39  displays  [H()3)]”  against  ^  for  a  20-elcment  line  with  three  signals  in  O-dB 
uncorrelated  noise.  The  independent  variable,  as  in  the  case  of  the  circle  array,  is  the  array 
parameter  or,  equivalently,  the  Hue  length.  Figures  40  through  43  are  representative  ME  array- 
response  plots  for  random  lines  of  various  lengths.  The  signal  plus  noise  field  consists  of  a 
20-dB  SNR  signal  at  0  =  0°,  a  10-dB  SNR  signal  at  0  =  45°,  and  a  0-dB  SNR  signal  at  0  =  90°. 
Figure  40  also  displays  the  MUSIC  beamformer  and  conventional  array  response. 
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VIII.  DETECTION  PROBABILITIES  FOR  MINIMUM-ENERGY 
PROCESSING  WITH  RANDOM  ARRAYS 


The  peak  spurious  peak  impacts  the  probability  of  detecting  a  weak  signal  in  an 
obvious  way.  In  this  section,  we  evaluate  the  probability  of  detecting  a  signal  when  minimum 
energy  processing  is  used  with  a  random  array. 

Suppose  that  the  noise  field  is  described  by  a  correlation  matrix 


R  =  N(^  SqS*  +  Q 


(34) 


where  S^  is  the  signal  we  are  interested  in  detecting  and  Q  describes  the  remainder  of  the  noise 
field.  In  particular,  suppose  the  background  noise  field  consists  of  K  interfering  signals  S[,. . . , 
Sj(  of  powers  aj,. respectively  in  uncorrelated  noise  of  power  <7^,  Then, 


Q  =  No^  SjSj  +  ...  +  S|^SJ  +  oj  j 


(35) 


Here  we  assume  that  the  interest  signal  Sg  is  in  the  region  away  from  the  K  interferers. 
With  R  as  described  by  (34),  a  useful  fact  is  that 


R~*  = 


Q-l  Sg  Sg  Q-‘ 
S^Q-^Sg  +  I/No2 


The  minimum-energy  filter  vector  Wg  when  the  array  is  stored  in  the  direction  of  the  signal  Sg 
is  then 


R~*  S, 


Wo- 


S^R-lSg 


Q-‘  SgS*Q-‘ 
sjQ-'Sg  +  i/No| 


1-1 


Q->  SgS^Q  * 


SoQ'‘s„. 


'iH 


Q-«  S, 


SgS;Q-‘s„ 


sJq-'s 


s’Q-^SgSjQ-^S, 


sjQ-'Sg  +  l/No| 
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Q  'So 
s^q'So 

Q  '  So 


s^o's,, 

sJq-'Sq  +  l/No| 

sJq'So 

SoQ  '  So  +  l/Nol 


SoQ-'s, 


Hence,  the  output  power  of  the  ME  beamformer  in  the  direction  of  the  signal  of  interest  is 

Pme  =  <  Rw. 


So  Q’* 


^oV  -  r  ,  1  So 

— _  N,^  S„  S„  .  Q 

oO’  So  L  J  SoQ-'So 


NoJ  + 


(36) 


SoQ'‘So 


Since  Q  contains  no  signals  in  the  direction  of  So,  our  earlier  sidelobe  approximation  is 
valid  for  the  second  term  in  (36).  Writing 

1 

SoQ^'So 


we  have  that 


K 

1  -  X  IS^Sol^ 

L=1 


If  the  array  is  random,  y  is  a  random  variable  with  a  probability  density  function 


h(y)  =  - 

y2(k  -  1)!((t2)‘^ 


(37) 
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This  is  just  equation  (28)  with  y  in  place  of  p.  The  output  power  p  in  the  direction  of  the  inter- 
est  signal  Sq  is  then  Na^  +  y  and  is  a  random  variable  with  a  probability  density  function 


j(P) 


(p  -  N<72)2(K  -  l)!(a2)K 


p  -  Na|  _j 


K-1 


exp 


1  - 


p  -  No|  ^ 


’  P  -  ^ 


(38) 


There  are  K  +  1  signals  in  the  noise  field  counting  the  interferers  and  the  signal  of  inter¬ 
est.  In  the  region  away  from  any  of  these  signals,  the  array-response  level  is  a  random  variable 
with  a  probability  density  function  similar  to  h(y)  only  with  K  +  I  replacing  K.  It  is 


e(y) 


4 


y2K!(<72)K+l  1_ 

and  the  corresponding  distribution  function  is 


■  ^  ■ 

K 

1 

V 

1 

y 

exp 

/  7 

1 

■  JJ 

•  V  >  ^2 

.  y  —  t'n  ’ 


Uy)  =  f  e(t)dt 


1  -  exp 


4  nr, 

K!  L„2  l 

yJJ  J 

(39) 


Compare  this  with  equation  (28). 

As  presented  in  the  previous  section,  the  peak  spurious  peak  is  the  largest  of  the  n  inde¬ 
pendent  samples  of  the  array-response  pattern  where  n  is  the  array  parameter.  It  is  a  random 
variable.  Let  r  represent  the  level  of  the  peak  spurious  peak  and  m(r)  be  its  probability  density 
function.  The  density  m(r)  is  unknown,  but  the  distribution  function  M(r)  is  found  to  be 

r  n 

M(r)  =  f  m(t)dt  =  [  Ur)  ]  .  (40) 
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See  equation  (31). 

For  probability  of  signal  detection,  we  take  the  probability  that  the  signal  peak  p  will 
e.Kceed  the  peak  spurious  peak  r.  Assuming  that  p  and  r  are  independent,  this  is 


Pr[p  >  r]  =  j(p)n(p)drdp 

p>r 


0) 


=  /  [  f  m(r)dr  ]  j(p)dp 


No2  +  a2  aj 


OJ 

/ 

Na2  +  al 


M(p)j(p)dp 


(O 


-  j  CMp)]  j(p)dp 


(41) 


No|  +  oJ 

where  j(p)  is  given  by  (38)  and  L(p)  by  (39). 

Making  the  substitution  =  l/N  and  the  change  of  variables  x  =  l/o^  (1  -  o^/p), 
this  integral  becomes 


w 


Pr[p  >  r]  =  f  [  L(p)  ]  j(p)dp 


N(7|  +  al 


N 


=  / 


KtIt) 


1  -  e 


-X 


N" 


(N  -  5(N  -  x))2(K  -  I)! 
r  /(6  +  l)x  -  N6\  1 ) 

*  L  ^  \5x  +  (I  -  5)n/  J  ) 


1  +  X  +  +  — 

k! 


r  (0 .  ■„  -  N.\-| 

,  [  \  fix  Ml  -  fi)N/_ 


K-1 


dx 
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with  6  =  a^.  This  integral  was  evaluated  numerically. 

Figure  44  is  a  plot  of  this  probability  against  input  SNR  (Og/ a^)  for  N  =  20,  n  =  5 1, 
and  K  =  2.  This  corresponds  to  signal  detection  with  a  20-element  random  circular  array  of 
diameter  D  =  10  X  in  the  presence  of  two  interferers. 

Following  figure  44  is  a  series  of  plots  displaying  precisely  this  situation.  A  series  of 
20-element  random  circular  planar  arrays  were  generated.  The  “background”  noise  field 
consists  of  0-dB  white  noise  and  two  interferers.  The  interferers  are  located  at  0  =  45°  and  90° 
and  each  has  a  SNR  of  3  dB.  The  signal  of  interest  —  the  one  we  want  to  detect  —  is  located 
at0=  180°. 

In  figure  45,  the  signal  of  interest  has  a  SNR  of  -10  dB.  Referring  to  figure  44,  we  see 
that,  at  this  SNR,  there  is  an  extremely  high  probability  that  the  signal  peak  will  exceed  the 
peak  spurious  peak.  Indeed,  the  signal  of  interest  is  easily  detected  using  minimum  energy 
processing. 

In  figures  46  through  48,  the  interest  signal  has  a  SNR  of  -13.2  dB.  At  this  source  level, 
the  probability  is  —  .9  that  the  signal  peak  will  exceed  the  peak  spurious  peak.  For  these  three 
random  arrays  at  this  source  level,  the  signal  is  easily  detected  as  we  would  expect. 

In  figures  49  through  51,  the  interest  signal  has  a  SNR  of  -15.5  which  corresponds  to 
a  probability  of  —  .5.  Here,  we  see  that  the  signal  is  being  lost  in  the  side  region  noise.  In 
figure  49,  the  signal  peak  is  about  equal  in  height  to  the  peak  spurious  peak.  In  figure  50,  for 
another  realization  of  the  random  array,  the  signal  peak  is  the  highest  peak.  In  figure  51,  the 
signal  is  clearly  lost  in  the  noise. 

The  final  figure  of  this  series,  figure  52,  shows  the  minimum-energy  array  response  for 
this  noise  field  and  a  representative  array  when  the  SNR  of  the  interest  signal  is  -18.0  dB.  This 
corresponds  to  a  probability  of  =  .  1  and,  as  expected,  the  signal  is  lost  in  the  noise. 

Figures  53  and  54  are  similar  to  figure  44.  Both  show  the  probability  that  the  signal 
peak  exceeds  the  peak  spurious  peak  (i.e.,  the  detection  probability)  for  20-element  random 
arrays  against  input  SNR.  The  first,  figure  53,  is  for  a  noise  field  containing  two  interferers. 
Curves  (left  to  right)  are  for  n  =  31, 41,  51,  71,  82,  and  103  corresponding  to  circular  planar 
arrays  of  diameters  D  =  6X,  8X,  lOX,  14X,  16X,  and  20X  respectively.  Note  that  for  a  fixed 
number  of  elements,  increasing  the  physical  size  of  the  array  decreases  the  probability  of  signal 
detection. 

Figure  54  is  for  a  noise  field  containing  three  interferes.  Curves  (left  to  right)  corres¬ 
pond  to  array  diameters  of  6X,  8X,  lOX,  14X,  and  20X.  Comparing  figures  53  and  54,  we  see 
that  increasing  the  number  of  interferers  complicates  the  noise  field  and,  consequently,  the 
array-response  pattern  making  detection  of  weak  signals  more  difficult. 
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Figure  20.  AAR  and  CBP.  20-element  random  circle;  diameter  =  6  wavelengths. 
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Figure  22.  Deterministic  positions.  20-element 
rectangular  array. 


Figure  23.  Conventional  beampattern.  20-element  rectangular  array. 


Figure  24.  ME  array  response.  Rectangular  array;  one  signal. 
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Figure  32.  Peak  spurious  peak  estimate.  Random  circle  diameter  =  10  wavelengths; 
N  =  20  K  =  3. 
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Figure  34.  Random  circular  planar  array.  20-elements;  various  diameters. 
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Figure  38.  ME  array  response.  Random  circular  array  diameter  =  14  wavelengths. 
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Figure  40.  ME  array  response.  20-etement  random  line;  length  =  10  wavelengths. 
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Figure  41.  ME  array  response.  20-etement  random  line;  length  =  20  wavelengths. 
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Figure  42.  ME  array  response.  20-element  random  line;  length  =  30  wavelengths. 
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Figure  45.  ME  array  response.  Interest-signal  SNR  = -10.0  dB. 
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Figure  46.  ME  array  response.  Interest-signal  SNR  =  -13.2  dB,  case  1. 
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Figure  53.  Detection  probability.  Random  circular  array;  N  =  20,  K  =  2. 
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Figure  54.  Detection  probability.  Random  circular  array;  N  =  20,  K  =  3. 
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APPENDIX  I 


Theorem 


If  1|P||  <  l,then(I-P)'*  exists  and  ([I-O-P)  *11  -  — 


IIPI! 


Proof: 


First,  let  x  0.  Then 

ll(I-P)llx||  =  ||x-Px||>||x|t-||Pxli>||x||-||P||  .  ||x|| 

=  (1  -  IIPID-  llxl!>0 
so  that  (I  -  P)x  #  0  and  (I  -  Pr  *  exists. 

1 


Next,  we  will  show  that  |  |(I  -  P)  M  I— 


1-tlPl 


From  (I  -  P)  (I  -  P)- '  =  (I  -  P)"*  -  P(I  -  P)  *  =  I  . 
we  have  that  (I  -  P"^  =  I  P(I  -  P)”* 

.•.ii(i-  p)-Mi  <  mil  +  iiPii  •  11(1-  p)-*  II  = » ^  iipii  •  IK*  - 
(i-iiPii).ii(i-pr'(i^i 
1 


ll(i-pr‘ll^ 


i-llPlI 

Finally,  from  (I  -  P)  (I  -  P)"*  =  d  -  P)'*  -  Pd  -  P)-'  =  I  , 
we  have  that  I  -  (I  -  P)->  =  -  P(I  -  P)~* 

IIPII 


.•.lli-d'Pr*ll^llPll-ll(*-pHll^ 


I -IIPII 
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APPENDIX  II 


Theorem 

Let  X|X2,  . . X^  be  i.i.d.  random  variables  with  a  probability  density  function 
f(Xj)  =  —  e  '  Xj  >  0  ;  i  -  1,2,  K  . 

Then  Y  =  Xj  +  X2  +  . . .  +  Xj^  has  a  distribution  with  a  p.d.f. 

g(y)  =  - ;  y  >  0  . 

(K  -  1)!(o2)K 

Proof: 

The  Laplace  transform  of  one  of  the  Xj  s  is 


Lr(f(Xi)) 


00 

e  f(Xi)dxj 
^  0 
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The  Laplace  transform  of  X  is 


L^g(y))  =  L^f(Xi))  =  - -  =  [1  *  (A- 

L=1  L=1  1  +  o^r 

Here,  Lj(g(y))  is  analytic  in  the  strip  corresponding  to  R^r]  >  -l/o^.  Rather  than  take  the 
universe  Uplace  transform,  we  just  take  the  Uplacc  transform  of  the  p.d.f.  in  the  proof  state 

ment.  This  is 


L;(g(y))  = 


*  yK-lg-y/o^dy 
(K  -  l)!(a2)K 


(K  - 


Integrating  by  parts  with  u  = 


,K-I 


t-i] 


(K  - 


yielding 
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Continuing  in  this  fashion, 


t;(g(y))  = 


OO 


=  [1  .  ra^T^ 


Noting  the  equality  between  (A-1)  and  (A-2)  we  now  have  the  result. 


(A-2) 
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